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The 15 0(a,7) 19 Ne reaction is one of two routes for breakout from the hot CNO cycles into the rp 
process in accreting neutron stars. Its astrophysical rate depends critically on the decay properties 
of excited states in 19 Ne lying just above the 15 + a threshold. We have measured the a-decay 
branching ratios for these states using the p( 21 Ne,t) 19 Ne reaction at 43 MeV/u. Combining our 
measurements with previous determinations of the radiative widths of these states, we conclude 
that no significant breakout from the hot CNO cycle into the rp process in novae is possible via 
15 0(a,7) 19 Ne, assuming current models accurately represent their temperature and density condi- 
tions. 

PACS numbers: 26.30.+k, 25.60.Je, 26.50.+X, 27.20. +n 



Novae are thermonuclear runaways initiated by the ac- 
cretion of hydrogen- and helium-rich material from stellar 
companions onto the surfaces of white dwarfs in binary 
systems. Energy production and nucleosynthesis in the 
hottest novae are determined principally by the CNO, 
NeNa, and MgAl cycles jj] . Under high temperature and 
density conditions, e.g. in accreting neutron stars, break- 
out from the hot CNO cycles into the rp process occurs 
dramatically increasing the luminosity of outbursts 
and synthesizing nuclei up to masses of 100 u j|. Sev- 
eral reactions have been suggested as pathways for this 
breakout but only two are currently thought to be 
possibilities: 15 0(a,7) 19 Ne and 18 Ne(a,p) 21 Na. In as- 
trophysical environments the 15 0(a,7) 19 Ne reaction pro- 
ceeds predominantly through resonances lying just above 
the 15 + a threshold at 3.529 MeV in 19 Ne, as the 
direct capture component is very small by comparison 
|^|, p|. The reaction rate in novae is determined by the 
a-width T Q of the 4.033 MeV 3/2+ state, owing both to 
its close proximity to the 15 O + a threshold and its low 
centrifugal barrier to a-capture. 

A previous attempt to determine T a for this state was 
based upon measurements of a-transfer reactions to the 
analog state in the mirror nucleus 19 F j7). Such deter- 
minations, however, are subject to large uncertainties 
PL Direct measurements of the low energy cross sec- 
tion, which require high- intensity radioactive 15 O beams, 
are planned. At 19 Ne excitation energies relevant to no- 
vae and accreting neutron stars, only the a- and 7-decay 
channels are open, as the proton and neutron separation 
energies are 6.4 and 11.6 MeV H respectively. Hence, 
by populating these states and observing the subsequent 
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a- and 7-dccays, one can deduce the branching ratio 
B Q = r Q /r. If T 7 is also known, one can then calcu- 
late r a and thereby the contribution of each state to the 
resonant rate of 15 0(cv,7) 19 Ne. A pioneering effort of this 
kind was made by detecting a particles from the decay 
of 19 Ne states populated via the 19 F( 3 He,i) 19 Ne reaction 
JTo| | , but the sensitivity of the experiment was insufficient 
to measure B Q for the critical 4.033 MeV state, which 
was expected to be of order 10~ 4 . Despite vigorous ef- 
forts worldwide ]Tl] , p"2| , up to now no experiment has 
reached this level. 

In an experiment at the Kernfysisch Versneller Insti- 
tuut, we have obtained branching ratio data of high sensi- 
tivity by applying a novel method introduced at Argonne 
National Laboratory using a different reaction |0 . Pop- 
ulating the important states via the 21 Ne(p, ip^Ne re- 
action in inverse kinematics with a 21 Ne beam energy 
of 43 MeV/u, we detected either 19 Ne recoils or their 
15 O a-decay products in coincidence with tritons in the 
Big-Bite Spectrometer (BBS) [O. The large momentum 
acceptance of the BBS (Ap/p = 19%) allowed detection 
of either 19 Ne recoils or 15 decay products along with 
tritons emitted backward in the center of mass system. 
Positioning the BBS at 0° maximized the yield to the 
4.033 MeV 3/2+ state in 19 Ne. This state, whose dom- 
inant shell-model configuration is (sd) 5 (lp)~ 2 was 
selectively populated by an £ = 0, two-neutron transfer 
from the 3/2 + ground state of 21 Ne. Position measure- 
ments in two vertical drift chambers (VDCs) [|l6| allowed 
reconstruction of the triton trajectories. Excitation ener- 
gies of the 19 Ne residues were determined from the kinetic 
energies and scattering angles of the triton ejectiles. The 
7-decays of states in 19 Ne were observed as 19 Ne-triton 
coincidences in the BBS, whereas cv-decays were identi- 
fied from 15 0-triton coincidences. 

Recoils and decay products were detected and stopped 
just in front of the VDCs by fast-plastic/slow-plastic 
phoswich detectors [p"7| that provided energy loss and 
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FIG. 1: 19 Ne-triton coincidences (7-decays of states in 19 Ne). 
The curve is the sum of a constant background and 8 Gaus- 
sians centered at the energies of known states in 19 Ne ^j, the 
widths of which were determined by the experimental resolu- 
tion of 90 keV FWHM. 
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FIG. 2: 15 0-triton coincidences (a-decays of states in 19 Ne). 
The curve is the sum of a constant background and 6 Gaus- 
sians corresponding to known states in 19 Ne, the widths of 
which were fixed by the experimental resolution. The 15 + 
a threshold lies at 3.529 MeV. 



total energy information. A separate array of phoswich 
detectors was used to identify tritons after they passed 
through the VDCs. Timing relative to the cyclotron ra- 
dio frequency signal was also employed for unambiguous 
particle identification. An Al plate prevented many of the 
heavy ions copiously produced by projectile fragmenta- 
tion reactions of the 21 Ne beam in the (CH2)„ target from 
reaching the VDCs. The spatial extent of the heavy-ion 
phoswich array was sufficient to guarantee 100% geomet- 
ric efficiency for detection of 19 Ne recoils and 15 O decay 
products for 19 Ne excitation energies < 5.5 MeV. This 
resulted largely from the forward focusing of the 19 Ne 
recoils, which emerged at angles < 0.36° for tritons with 
scattering angles < 4°. The low decay energies of the 
states studied limited the angular and energy spreads of 
the 15 O decay products. High geometric efficiency, com- 
bined with the excellent background rejection afforded 
by detecting 15 nuclei instead of a particles, provided 
sensitivity to very small a-decay branches. 

The 19 Ne excitation energy spectrum obtained from 
19 Ne-triton coincidences, representing 7-decays of states 
in 19 Ne, is shown in Fig. |l|. Its most prominent peak 
is due to the 4.033 MeV 3/2+ state. The 4.379 MeV 
7/2+ state is the second most strongl y p opulated. Con- 
tributions from other known states |18| are indicated. 
The experimental resolution of 90 keV FWHM is insuf- 
ficient to resolve the 4.140 and 4.197 MeV states from 
one another; the 4.549 and 4.600 MeV states are also un- 
resolved. However, the astrophysically important 4.033 
and 4.379 MeV states are well separated from the oth- 
ers. The curve shown in Fig. [l] is the sum of a constant 
background and 8 Gaussians centered at the known ener- 
gies of the states, with standard deviations fixed by the 
experimental resolution of 90 keV FWHM. 

Fig. [2] shows the 19 Ne excitation energy spectrum ob- 
tained from 15 0-triton coincidences, corresponding to a- 
decays of states in 19 Ne. The observed states are labeled 



by their energies; the 4.549 and 4.600 MeV states can- 
not be resolved. A fit consisting of Gaussians plus a 
constant background is shown as well. The states below 
4.549 MeV decay overwhelmingly by 7-emission, while 
the higher-lying states observed here decay preferentially 
by a-emission, simply because the larger available de- 
cay energy enhances the barrier penetrability. The back- 
ground represents a larger fraction of the total events in 
the 15 0-triton coincidence spectrum than in the 19 Ne- 
triton coincidence spectrum, though it is still very small. 
This background is well reproduced by a constant, the 
value of which is determined to a precision of 7% (ler). 

No statistically significant evidence for a-decays from 
the states at 4.033 and 4.379 MeV was observed. For 
these states the a- and 7-decay spectra were numerically 
integrated in 100 keV intervals centered at the known en- 
ergies of the states. These data were subjected to both 
Bayesian and classical statistical analyses to determine 
upper limits on the a-decay branching ratios at various 
confidence levels. The two analyses agreed rather well, 
with the calculated upper limits differing by less than 
20% in all cases. The Bayesian analysis was found to 
be more conservative, and has been adopted here. As 
expected, there is no indication in the data of a-decays 
from the 4.140 and 4.197 MeV states because these de- 
cays are hindered by £ = 4 centrifugal barriers and low 
decay energies. For the states at 4.549, 4.600, 4.712, 
and 5.092 MeV, B a was obtained from the fits described 
above. The branching ratios are shown in Table | along 
with the results of Refs. Jl^, [l{J. Uncertainties in the 
present branching ratio determinations are purely statis- 
tical. Our 90% confidence level upper limit on Bq, for the 
4.379 MeV state is a factor of 11 smaller than the central 
value of Ref. (To) , a discrepancy we attribute to imperfect 
background subtraction in the previous determination. 
To obtain the resonance strengths we take a weighted 
average of B Q where more than one measurement is avail- 
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able, excepting the 4.379 MeV state for which our upper 
limit is preferred. The adopted values are shown in Table 
H The uncertainties given in the table are la values, and 
we specify all upper limits at the 90% confidence level. 

The experimental data on T 7 for states in 19 Ne are 
sparse. Of the six states considered here, measurements 
are available only for the 4.033 MeV state. The value 
adopted for this state |2Q] is the result of a combined 
analysis of Coulomb excitation and Doppler shift atten- 
uation J2][ data using shell-model calculations of the rel- 
ative strengths of E2 and Ml transitions. In some cases, 
widths of the analog state from the mirror nucleus 19 F 
have been measured, and we adopt these under the as- 
sumption that r 7 ( 19 Ne) = T 7 ( 19 F). Such measurements 
are available for the 4.549 [[§], 4.600 |22j, and 4.712 MeV 
states @. For the 4.379 and 5.092 MeV states, measure- 
ments in neither nucleus are available, and we adopt the 
results of shell-model calculations |25|] , assigning a la un- 
certainty of 20% to the calculated widths. The values of 
r 7 and r Q , which is calculated as T a — 1 B g T 7 , are 
shown in Table |. For the 4.033 and 4.379 MeV states we 
calculate upper limits on r a at the 90% confidence level 
using 1.64cr upper limits on T 7 . Earlier compilations of 
decay widths can be found in Refs. || p4| . 
We have calculated 
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thermally averaged 
reaction rates due to these 6 resonances. Both the in- 
dividual rates and the sum of the resonant and direct 
capture contributions are shown in Fig. |3[ The contri- 
butions of the 4.033 and 4.379 MeV states are calculated 
using our 99.73% confidence level upper limits for their 
a- widths, 31 /ieV and 5.6 meV respectively. We do not 
show the individual contribution of the 4.549 MeV state 
because it is insignificant by comparison with the other 
resonances. The direct capture rate was calculated as in 
Rcf. [§), but is significant only below 0.1 GK. Our upper 
limit on the contribution of the 4.033 MeV state is much 
larger than all other contributions to the reaction rate 
for T < 0.5 GK. On the contrary, the 4.600 and 4.712 
MeV states account for most of the reaction rate at the 
high temperatures of 1.9 GK found in accreting neutron 
stars ||. 

The amount of leakage from the hot CNO cycle via 
15 0(a,7) 19 Ne depends on its rate compared to the /3 + 
decay rate of 15 (ti/2 = 122 s). In order to calculate 
the reaction rate in a particular environment one needs 
to know the local He mass fraction Y. We assume here a 
Y of 0.27, which is the solar value |26| adopted in accret- 
ing neutron star models but is approximately twice 
the maximum value used in nova models (2?], g9| . Al- 
though the accreted material in novae is usually assumed 
to be of solar composition, significant mixing with the 
surface material of the white dwarf occurs prior to the 
nova outburst |3(J , rendering the net Y smaller than that 
of the accreted matter. For this reason and because we 
adopt a 99.73% confidence level upper limit for the con- 
tribution of the 4.033 MeV state, the reaction rate we 
calculate with this Y represents an extreme upper limit 
for novae. Fig. |^ shows the boundary in the density- 




FIG. 3: Product of the Avogadro constant Na and the ther- 
mally averaged rate of the 15 0(a,7) 19 Ne reaction per particle 
pair. Contributions from all of the important resonances are 
shown, along with the sum of the resonant and direct capture 
rates. The contributions of the 4.033 and 4.379 MeV states 
are 99.73% confidence level upper limits. 



temperature plane at which the rate of the 15 0(a,7) 19 Ne 
reaction equals the (3 + decay rate of 15 for a Y of 0.27. 
Also shown are shaded regions corresponding to the peak 
temperatures and densities reached in nova outbursts and 
accreting neutron stars. A comparable figure appears in 
Ref. 

Typical nova models reach peak temperatures of 0.2 - 
0.3 GK, at which time densities are of order 100 g cm~ 3 
pSl p9 . 32 1 . Under such conditions, the /3 + decay rate of 
is more than 3 orders of magnitude faster than the 
rate of the 15 0(a,7) 19 Ne reaction. At temperatures be- 
low 0.1 GK in novae, energy is generated by the pp chains 
and the cold CNO cycle. Only at temperatures above 
0.1 GK can the 13 N(p, 7) 14 reaction compete equally 
with 13 N (3 + decay, initiating the hot CNO cycle. Dur- 
ing a nova outburst the temperature exceeds 0.1 GK for 
about 1000 s ||^, |32|) , while the mean time required to 
complete one loop of the hot CNO cycle, determined by 
the /?+ halflives of 14 and 15 and the 13 N(p, 7 ) 14 
rate, is at least 300 s. Hence no more than a few cycles 
can be completed during a nova explosion, and the frac- 
tional leakage per cycle through 15 0(a,7) 19 Ne of < 0.001 
cannot process a significant amount of CNO material to 
higher masses. Since the rate of the 18 Ne(a,p) 21 Na re- 
action appears far too small in novae to compete with 
18 Ne j3 + decay j|, |H, [?4|, we conclude that appreciable 
breakout from the hot CNO cycle into the rp process in 
novae is precluded given our current knowledge of reac- 
tion rates and nova physics. This conclusion is consistent 
with those reached in a recent study of reaction rate vari- 
ations on nova nucleosynthesis p5[ , which considered a 
range of 15 0(a, 7 ) 19 Ne rates from 0.002 to 30 times the 
rate used here. 



In summary, we have measured the a-decay branch- 
ing ratios for all of the states in 19 Ne relevant to the 
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TABLE I: Branching Ratios B a = T a /r and Decay Widths. Upper limits are specified at the 90% confidence level. 



E x (MeV) 



B a (present work) B a (Ref. 1C]) B Q (Ref. [19]) B a (adopted) P 7 (meV) Ref. 
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< 4.3 xl0~ 4 

< 3.9 xlO -3 
0.16 ± 0.04 
0.32 ± 0.04 
0.85 ± 0.04 
0.90 ± 0.06 



0.044 ± 0.032 
0.07 ± 0.03 
0.25 ± 0.04 
0.82 ± 0.15 
0.90 ± 0.09 



0.32 ± 0.03 



< 4.3 xlO" 4 

< 3.9 xl0~ 3 
0.10 ± 0.02 
0.30 ± 0.02 
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196 ± 39 
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FIG. 4: Density at which the rate of the 15 0(a,7) 19 Ne reac- 
tion equals the (3 + decay rate of 15 O for a He mass fraction 
of 0.27, the solar value. Below 0.5 GK, the curve represents 
our 99.73% confidence level lower limit. The shaded regions 
indicate the peak temperature and density conditions found 
in novae and accreting neutron stars. 



astrophysical rate of the 15 0(a,7) 19 Ne reaction, popu- 



lating them by means of the p( 21 Ne,i) 19 Ne reaction at 43 
MeV/u and observing their decays with 100% geomet- 
ric efficiency in a magnetic spectrometer. Combining our 
measurements with prior determinations of the 7-widths 
of these states, we have calculated the astrophysical rate 
of 15 0(a,7) 19 Ne. For the first time, the calculation of the 
reaction rate at T < 0.5 GK is based on direct experimen- 
tal measurements. On the basis of these calculations and 
independent determinations of the 18 Ne(a,p) 21 Na reac- 
tion rate, we conclude that there can be no appreciable 
breakout from the hot CNO cycle into the rp process 
in novae, assuming current models accurately represent 
their temperature and density conditions. 
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